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ABSTRACT 

If cluster dark matter is in the form of compact objects it will introduce fluctuations into the 
light curves of distant sources. Current searches for MACHOs in clusters of galaxies focus on 
monitoring quasars behind nearby systems. This paper considers the effect of such a compact 
population on the surface brightness distribution of giant gravitationally lensed arcs. As the 
microlensing optical depth is significant in these clusters, the expected fluctuations are substantial 
and are observable. Focusing on the giant arc seen in Abell 370, we demonstrate that several 
'extreme' events would be visible in a comparison of HST observations at two epochs. Utilizing 
NGST, long term monitoring should reveal a ubiquitous twinkling of brightness over the surface 
of the arcs. 

Subject headings: galaxies: clusters: general - gravitational lensing - dark matter 



1. Introduction 

After many years of searching, the nature of 
dark matter still eludes us. Over the last decade, 
however, programs searching for microlensing vari- 
ability in stars in the Magellanic Clouds have re- 
vealed a population of compact objects in front 
of the LMC (Afonso et al. 1999; Alcock et al. 
2000a). The EROS team have used their obser- 
vations to place an upper limit of 10% of the 
mass fraction in MACHOs in the Galactic halo 
(Lasserre et al. 2000), while the MACHO team 
interprets their data as providing evidence of a sig- 
nificant fraction, perhaps 20%. The recent discov- 
ery of nearby cold, high proper motion, hydrogen 
atmosphere white dwarfs appears to support the 



latter interpretation (Ibata et al. 2000; Hodgkin 
et al. 2000). These results are complicated by 
fact that views towards the Galactic Bulge ap- 
pear to be overdense in compact objects (Alcock 
et al. 2000b), suggesting that the Galactic Halo 
is clumped (c.f. Klypin et al. 1999) and our view 
towards the Magellanic Clouds simply represents 
an underdense line-of-sight. Currently, therefore, 
the true fraction of the Galactic Halo composed of 
compact objects is uncertain. 

Galaxy clusters represent the largest bound ac- 
cumulations of matter in the Universe, and a 
number of lines of evidence, including kinematic 
(Bird, Dickey & Salpeter 1993), X-ray (Henry et 
al. 1993) and gravitational lensing (Kneib et al. 



1993) studies, indicate that they must possess a 
substantial quantity of dark matter. In this paper 
we address the question of what would be the ob- 
servational consequences if this dark matter were 
also composed of compact objects. 

An approach to test the nature of cluster dark 
matter was recently suggested by Walker & Ire- 
land (1995) and Tadros, Warren k Hewett (1998; 
2000). As with compact matter in the Galac- 
tic Halo, MACHOs in clusters are expected to 
introduce brightness fluctuations in background 
sources, and this is searched for in a monitor- 
ing program of ^ 600 quasars behind the Virgo 
cluster. While initial results are promising, the 
vicinity of the Virgo cluster means that the mi- 
crolensing optical depth is small 0.001) and in- 
duced fluctuations are relatively rare. The optical 
depth can be increased by changing the lensing 
geometry and considering clusters at higher red- 
shift. With this, however, it is difflcult to identify 
a substantial number of background quasars with- 
out greatly increasing the number of clusters to 
be monitored. The problem is also compounded 
by the difficulty of distinguishing microlensing- 
induced quasar variability from intrinsic mecha- 
nisms. 

Lewis & Ibata (2000) addressed the question 
of what efii'ect a cosmological distribution of com- 
pact objects would have on the surface brightness 
distributions of galaxies at z < 0.5. While the mi- 
crolensing optical depth is quite small (cr ^ 0.04), 
the resulting low-level fluctuations, ^ 2%, are ob- 
servable. In this paper we extend this analysis 
to the view of distant galaxies observed through 
galaxy clusters whose dark matter is composed of 
MACHOs. The ideal targets are the giant gravi- 
tationally lensed arcs. Focusing on these systems 
ensures that the column density of matter through 
the cluster is substantial, resulting in an optical 
depth near unity; this vastly exceeds 'self-lensing', 
with an optical depth of only 10^^, where MA- 
CHOs microlens galaxies within the cluster (Gould 
1995). Potential sources abound, as giant arcs 
and extended images have been identified in more 
than thirty clusters (e.g Fort & Mellier 1994). As 
the induced variability will be searched for on a 
pixel-by-pixel basis, each extended image presents 
a large number of resolution elements which can 
be viewed as separate sources, far outweighing the 
potential number of background quasars. 



In Section 2 we outline the numerical approach 
we adopt to tackle this problem, focusing on the 
simulation of microlensing light curves. As an 
illustration of the efficiency of microlensing 
mechanism for introducing variability, we present 
in section 3 a case study of the giant arc seen in the 
cluster Abell 370. This section also discusses the 
observational considerations for a monitoring pro- 
gram of giant arcs. The conclusions to this study 
are presented in Section 4. 

2. Method 

In examining microlensing of sources in the 
Magellanic Clouds, the goal is to identify the vari- 
ation in brightness of a single star as a compact 
object passes in front of it. Considering instead 
microlensing towards more distant sources, indi- 
vidual stars can no longer be resolved and in- 
duced variations are identified against an unre- 
solved stellar background. Termed 'pixel-lensing' 
(e.g. Crotts 1992), the magnification of a star re- 
sults in the increase of fiux in an individual reso- 
lution element, and hence appears as a fiuctuation 
in the surface brightness of a source. Employing 
an efficient technique called 'Difference Imaging', 
searches for such fiuctuations have recently proved 
fruitful with the identification of microlensing vari- 
ability towards the Galactic Bulge and M31 (Al- 
cock et al. 1999a, 1999b; Ansari et al. 1999). 

In a similar vein, observations of gravitation- 
ally lensed giant arcs do not resolve individual 
stars, rather each pixel over an arc represents the 
summed flux from a population of stars. As stars 
in this population are magnified due to gravita- 
tional microlensing by compact objects in the fore- 
ground cluster, the apparent surface brightness of 
the pixels will be seen to fiuctuate. It is important 
to note that at the substantial optical depths con- 
sidered in this paper, each star in the population 
will be subject to microlensing and so the overall 
light curve seen will be the sum of the individual 
fluctuations of each star. This differs substantially 
from pixel lensing in the local group, where the 
optical depth is small and at most only a small 
number of stars are significantly microlensed at 
any instant. 

In simulating the expected variability, the fol- 
lowing procedure is adopted; firstly, for a partic- 
ular set of microlensing parameters, a large cat- 
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aloguc of light curves, fii{t), arc generated. A 
population of stars, Si, is randomly selected such 
that the total luminosity of the population is 
Lpop = S L{Si). Each star is used to normalize a 
light curve selected from the catalogue to produce 
a microlensed view of the star. These individual 
light curves are then combined to produce a to- 
tal microlensed light curve for the entire stellar 
population, M{t) = T, ^i{t) L{Si) . As a consis- 
tency check, the mean luminosity of the resultant 
light curves must equal {fJ'th)Lpop, where (fJ^th) is 
the mean magnification of the arc induced by the 
macrolcnsing model of the cluster. Typically, stel- 
lar luminosity functions increase at fainter magni- 
tudes, implying that populations are dominated, 
in number, by low luminosity stars. Such faint 
stars, even if substantially lensed, will not signifi- 
cantly contribute to the final light curve. Hence, 
for any star for which ^J.maxL{Si) / {mh)Lpop < 
0.005, where Umax is taken to be 70, the contri- 
bution to the final light curve is taken as uniform 
and equal to {iith)L{Si); such an approach reduces 
the required numerical calculations in construct- 
ing a light curve. As with Lewis & Ibata (2000), 
we adopt the stellar luminosity function as deter- 
mined by JahreiB & Wielen (1997). 

An efficient method of generating high optical 
depth microlensing light curves for a point-like 
source is employed (Witt 1993; Lewis, Miralda- 
Escude, Richardson & Wambsganss 1993). This 
approach treats the source trajectory as a line of 
infinite extent, ensuring that it is numerically sim- 
ple to find all the microlensed images with a con- 
tour following algorithm. This guarantees that 
the light curve is not lacking flux due to 'undis- 
covered' images. The efficiency allows large sam- 
ples of light curves to be generated, the basis for 
a number of extensive statistical studies (Lewis & 
Irwin 1995; Lewis & Irwin 1996; Wyithe & Web- 
ster 1999; Wyithe, Webster & Turner 1999,2000). 

While the approach is straight-forward, the vol- 
ume of potential parameter space is very large and 
a full examination of this is beyond the scope of 
this paper, and a more extensive study will be pre- 
sented elsewhere. As an illustration of the efficacy 
of microlensing on introducing variability into sur- 
face brightness distributions, we now focus on the 
spectacular giant arc in Abell 370. 



3. Case Study: A370 

At a rcdshift of z = 0.37, the rich galaxy cluster 
Abell 370 possesses several giant gravitationally 
lensed arcs (Soucail et al. 1987). Utilizing their 
spatial extent, detailed modeling of the mass dis- 
tribution of Abell 370 has been undertaken (Kneib 
et al 1993; Abdelsalam et al. 1998), confirming the 
binary nature of the cluster. The most prominent 
giant arc is over 20 arcsecs in extent and has a 
spectroscopic redshift of 2: = 0.724 (Soucail et al. 
1988). Examining the models for the mass distri- 
bution in Abell 370 in the vicinity of this giant arc, 
we choose a fiducial model with an optical depth 
of a = 0.6 and shear of 7 = 0.2. This results in a 
mean magnification of {fith) = 8.33. 

3.1. Light Curves 

A large sample of light curves was generated 
using the above parameter set and combined with 
the stellar luminosity function to produce the light 
curves of an entire population as outlined in Sec- 
tion 2. The source trajectory is aligned with the 
shear, although simple scaling relations to other 
shear orientations exist (Lewis & Irwin 1995). 
Each lens has a mass of IMq. 

Several light curves of these are presented in 
Figure 1 for stellar populations of Lpop = 10^, 
10'', 10^' and 10^L(5; these numbers represent the 
intrinsic luminosity of the population, the ob- 
served value is Lots = {^J■th)Lpop■ The abscissa 
is in units of Einstein radii for a solar mass star; 
this is related to a time scale, as explained below. 
The light curves display the complex variability of 
high optical depth microlensing. At low optical 
depths, the light curves would consist of solitary 
events due to individual microlenses. While these 
can be used to directly infer microlensing parame- 
ters, they occur very rarely and low optical depth 
light curves arc extremely quiescent. With the rich 
structure seen in the light curves of Figure 1, the 
details of the microlensing can be determined in a 
statistical sense (e.g. Lewis & Irwin 1996). 

In examining the light curves in Figure 1 it is 
apparent that the degree of variability drops as 
the total luminosity is increased, with populations 
of W^Lq displaying rapid variations of 10-20% on 
scales much less than an Einstein radius, while the 
light curves for populations of IO^Lq appear es- 
sentially flat with only small scale variability. This 
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is due to the fact that even a strongly lensed, very 
luminous star can contribute only a small amount 
to the total luminosity of a stellar population of 
Lpop ^ lO'^L©. 

The abscissa can be normalized into a time scale 
(Lewis & Ibata 2000); the crossing time of an Ein- 
stein radius of a Solar mass star is 9.2 I'looo'^"^ 
yrs (assuming Jlo = 1, Ao = 0), where wiooo is the 
velocity of the microlensing masses in units of 
1 000km/ s. When considering a population of mi- 
crolensing objects of a different mass it is only 
this crossing time that needs to be rescaled, with 
the degree of fluctuations in Figure 1 remaining 
the same. Being simply a function of the square 
root of the lensing masses, if the dark matter in 
Abell 370 was composed of Jupiter mass com- 
pact objects, the corresponding value would be 
~ 3.5W]^gQ/i^2 months. 

By combining the light curves from the various 
stellar samples it is possible to generate a prob- 
ability distribution for the expected fluctuations 
(Figure 2). These reflect the features observed in 
Figure 1, namely that microlensing of lower lu- 
minosity populations results in more substantial 
fluctuations. In terms of integrated probability, a 
lO^X© population will exhibit deviations from the 
mean exceeding 5% for half of all observations. 
A quarter of all observations will reveal fluctua- 
tions exceeding 10%. The corresponding values 
for the more luminous populations are 2% & 3.5% 
for Lpop = 10^, 0.6% & 1% for Lpop = 10*^ 
and 0.2% & 0.3% for Lpop = 10^. The Lpop = 
lO^L© case also displays more extreme fluctua- 
tions, experiencing deviations of > 15% for 10% 
of observations, and > 29% over 2% of its light 
curve. 

To detect the fluctuations expected for a popu- 
lation of cluster MACHOs it is advantageous to 
target stellar populations with smaller intrinsic 
luminosities. Targeting fainter sources, however, 
presents observational difficulties, such as substan- 
tially increased integration times, a point we turn 
to in the next section. 

3.2. Observational Considerations 

The study of cosmological microlensing by 
Lewis & Ibata (2000) examined stellar luminos- 
ity in a pixel covering a galaxy out to z ^ 0.5. 
Here we explore the possibility of using HST to 



constrain the microlensing variability in a strongly 
lensed arc of Abell 370. For a source redshift of 
z = 0.724, a spatially resolved population with in- 
trinsic U-band luminosity L = 10 ""Lq, magnified 
by a factor of 8.3, will give rise to an observed V- 
band flux of 1.3photons/sec/m^, (Oq = 1, Aq = 0, 
and Hq = lOOkm/s/Mpc). The expected total 
system efficiency of the new HST instrument the 
"Advanced Camera for Surveys" is approximately 

o 

40% at 6000 A, so 2.5photons/sec will be detected 

per IO^Lq population. The sky backgroimd will 
be negligible in comparison, O.lphotons/sec per 
resolution element (the area encompassed by the 
point spread function) ^. In a 10^ sec exposure 
HST will therefore measure variability to 0.66% 
accuracy; this would allow a 3a detection of 2% 
variability events which according to the simula- 
tions of §3.1, occur with 0.25% chance. The giant 
arc in Abell 370 subtends over ~ 1000 HST spa- 
tial resolution elements and 8 variability events 
would be detected in a 2-epoch observation cam- 
paign if all the cluster dark matter is in the form 
of MACHOs. Similarly, for a IQ^Lq population, 
the above 10'* sec integration time would allow 3% 
variability events to be detected at the 3(7 level. 
Approximately ~ 100 such events would be ex- 
pected in the 2-epoch dataset. Several exposures 
would be required at each epoch to correct for 
cosmic rays and to identify Poisson fiuctuations 
in the surface photometry that could be mistaken 
for microlensing events. 

Currently, only coarse surface photometry of 
the giant arc in Abell 370 has been published, but 
its total brightness suggests that each HST reso- 
lution element corresponds to an average intrinsic 
source luminosity of 10^ 10 ''L©. HST images 
of the arc reveals that it is very clumped (Small 
et al. 1996) and a range of source luminosities are 
available for study. We have begun a more de- 
tailed modeling of the microlensing in Abell 370, 
taking into account both the variation in the mi- 
crolensing parameters and surface brightness over 
the arc, with the goal of providing a map of the 
variability signature. 

To detect the more common variability in pop- 
ulations of lower luminosity will require NGST. 
With this instrument, the sky background is ex- 

^This value was derived from archival WFPC2 "Planetary 
Camera" data on this cluster. 
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pected to be substantially lower than on HST: < 
10~^ photons/sec/resolution element ^ at 1.1/xm. 
Assuming a 20% total system efficiency, a 10^ sec 
integration at two epochs would allow 10% vari- 
ability events to be detected at the 3a level in a 
population of intrinsic luminosity lO^i©. As dis- 
cussed above, variability of this magnitude must 
be very common, occurring in 25% of the resolu- 
tion elements per epoch. If MACHOs make up a 
substantial fraction of the dark matter in clusters, 
giant arcs will be seen to shimmer with NGST, 
the degree of shimmering being dependent upon 
the luminosity of the population under considera- 
tion. 

The above analysis is very encouraging, demon- 
strating that the signature of microlcnsing due to 
cluster dark matter could soon be detected with 
HST, as well as illustrating that microlensing fluc- 
tuations may be ubiquitous in future observations 
with NGST. 

4. Conclusions 

If dark matter in galaxy clusters exists in the 

form of compact objects, it should introduce vari- 
ability on background sources. Current searches 
for cluster dark matter have focused upon moni- 
toring quasars behind the Virgo cluster. Its close 
proximity, however, makes Virgo a poor lens. 

By considering clusters at higher redshift, the 
microlensing optical depth increases, although the 
potential number of background quasars that can 
be used as microlensing targets is reduced. In this 
paper wc have examined the microlcnsing effect 
that a population of compact dark matter objects 
would have on the surface brightness distributions 
of giant luminous arcs. Focusing upon the spec- 
tacular arc in Abell 370, numerical simulations re- 
veal that substantial variability can be introduced 
over the Icnscd image. This variability is, how- 
ever, dependent upon the intrinsic luminosity of 
the population of stars under consideration, with 
more luminous populations displaying more quies- 
cent light curves. 

To detect microlensing induced variability in 
the surface brightness of the giant arcs, a deep 
monitoring program over a single semester is re- 
quired. Such a program would allow easy identifi- 

^see http;//www.ngst. stsci.edu/sky/sky.html 



cation of variability contaminants, such as novae, 
supernovae and highly luminous variable stars, 
due to their characteristic temporal profiles (see 
Lewis & Ibata 2000). These contaminating sys- 
tems are also rare, while microlensing variability 
will appear ubiquitously over the arc. As the mass 
distribution in the cluster can be determined via 
the modeling of the macrolensing features, the 
identification of microlensing variability will re- 
veal the MACHO mass fraction within the clus- 
ter. Many strongly lensed systems are available 
for study, several displaying very extended struc- 
ture [e.g. cB58 (Yce et al. 1996)], providing ideal 
tools for the examination of the nature of cluster 
dark matter. 
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Fig. 1. — Two sets of examples of light curves for various stellar populations. The abscissa is in units of 
Einstein radii, while the ordinate presents the fluctuations with respect to the mean. Both sets display the 
same general trends that are seen in all the light curves. Namely, that the lower luminosity stellar populations 
show the most rapid variability. As the total luminosity of the population is increased, even luminous stars 
that are strongly lensed contribute only small variations to the light curve. 
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Fig. 2. — The fluctuation probability distributions for the populations with intrinsic luminosities of 
10^,10^,10^ and \^Lq. Note that the observed luminosity is Lohs = {lJ'th)Lpop and the different or- 
dinate scaling on each panel. As expected from the light curves presented in Figure 1, the lower luminosity 
populations show a broad distribution, which narrows as the total luminosity is increased. 
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